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Abstract 

Starting from a completely general standpoint, we find the most general 
brane-Universe solutions for a three-brane in a five dimensional spacetime. 
The brane can border regions of spacetime with or without a cosmological 
constant. Making no assumptions other than the usual cosmological symme- 
tries of the metric, we prove that the equations of motion form an integrable 
system, and find the exact solution. The cosmology is indeed a boundary 
of a (class II) Schwarzschild-AdS spacetime, or a Minkowski (class I) space- 
time. We analyse the various cosmological trajectories focusing particularly 
on those bordering vacuum spacetimes. We find, not surprisingly, that not 
all cosmologies are compatible with an asymptotically flat spacetime branch. 
We comment on the role of the radion in this picture. 
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I. INTRODUCTION 



The idea that our four-dimensional universe might in fact, upon closer inspection, be 
higher dimensional has always been compelling, although it has been conventionally sup- 
posed that the extra dimensions are small and compact. An almost diametrically opposed 
view, wherein the extra dimensions are large and noncompact has gained ground recently, 
in which our universe is a 'defect' (see for early work) in an anti-de Sitter bulk - the 
negative cosmological constant providing an effective or exotic compactification || of the 
spacetime. In most models, the universe is taken to be a domain wall, (i.e. of codimension 
one - although models of higher codimension have been considered [||) since this is most 
naturally supported by heterotic string theory ||-||. The original models of Randall and 
Sundrum |TJ focussed on the main features of having either one or two domain wall universes 
at orbifold fixed points of a single extra dimension; in the case of the two-wall model the 
second wall had negative energy and tension. The compelling feature of these brane-world 
models is that gravity has a four-dimensional character on the universe-domain-wall with 
short-range five-dimensional corrections coming from Kaluza-Klein (KK) modes. Indeed, it 
was demonstrated explicitly by Garriga and Tanaka || that linearized gravity on the brane 
was indeed Einstein gravity for a single wall universe (see also [|K||). The key feature of 



the Randall Sundrum (RS) (and indeed many other) models is that the induced metric on 
the wall (or walls) is flat Minkowski spacetime. In general, a domain wall spacetime is not 
static, but has a de-Sitter expansion along its spatial directions, however in the presence of 
a bulk cosmological constant, the effect of the wall energy and tension is neutralized, and 
a static solution is allowed. However, our universe is not Minkowski spacetime, nor is it 
a localized linearized perturbation thereof. Therefore, if brane universe models are to be 
relevant, cosmological brane solutions are the obvious next step. 

A general homogeneous and isotropic brane cosmology is simply a wall with energy, e, 
and tension, t, which are no longer fixed by the special domain wall relation e = t. For 
example, we might wish to set 

E = E + P ; T = T -p (1) 

and then take p and p to be our 'cosmological' energy density and pressure respectively. 
This is what is conventionally done, with p and p obeying an equation of state (rather than 
e and t), although it should be noted that when p and e are of a similar order of magnitude, 
there is no reason to suppose that the equation of state will continue to hold with this choice 
of p and p, which become rather arbitrary at that point. 

Much of the work on brane-cosmological models to date has taken a 'brane-based' ap- 



proach, notably the work of Binetruy et. al. |TT| , p!2[| (see also |L3|])) m which the RS spacetime 



is generalized to allow for time-dependent cosmological expansion. However, Binetruy et. al. 
were concerned primarily with deriving the four-dimensional brane-cosmological equations, 
and for simplicity when finding explicit solutions, chose a Gaussian Normal gauge in which 
g nn , the component of the metric normal to the wall, was set to unity. Therefore, while 
the brane cosmology was readily apparent in their approach, the bulk spacetime structure 
was less transparent, with coordinate singularities indicating the breakdown of the rather 
restrictive GN gauge. A more general brane-based approach is epitomized by the work of 



Shiromizu et. al. [TJj (see also Maartens, flT5|, for an excellent brane-based description of 
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bulk effects) in which the Gauss-Codazzi formalism is used to obtain a pure-brane 'Ein- 
stein' equation, in which the bulk geometry is encoded in a single tensor S^. While such 
a brane-based approach has the attraction of being a simple generalization of the stan- 
dard four-dimensional FRW equations, the effect of the bulk on the brane is somewhat less 
transparent, with the tensor £„„ hiding a multitude of sins. 



An alternate approach for deriving cosmological solutions was taken by Ida |16fl , who 
instead considered a 'bulk-based' point of view, in which the most general static AdS solution 
with the appropriate symmetries: 



ds 2 = h(r)dt 2 - h- l (r)dr 2 - r 2 



+ x 2 dn 2 H 



(2) 



(where h(r) = k — ^ + k 2 r 2 ) is taken, and the brane becomes an arbitrary boundary of 
two versions of this spacetime, with possibly differing masses on each side. The Israel- 
Gauss-Codazzi equations then give the energy and tension of the boundary as a function 
of its trajectory. The beauty of this approach is that it is quite general and does not 
rely on any Z 2 symmetry around the wall itself; it is also straightforward to derive brane- 
cosmological evolution equations - the cosmological solution then becomes a wall moving 
in AdS spacetime; the only disadvantage is that it is perhaps a little more abstract than 
the brane-based approach. A coordinate transformation relating the solutions of |ll],|12j to 
Ida's solutions was found by by Mukohyama et. al. [17]. Ida's work in fact generalizes the 



work of Kraus [TS[] (see also [Tj| for related string theoretic work, and [2D|, for an AdS/CFT 
perspective on the issue) , who derived the most general Z 2 symmetric domain wall solutions 
by taking slicings of (0), and so can apply to models with different cosmological constants 
such as the Lykken- Randall (LR) model, ||21|| , for example. 

What we aim to do in this paper is to bridge the gap between the 'brane-based' approach, 
where the brane represents a fixed boundary in some time-evolving spacetime, and the 'bulk- 
based' approach, where the wall (or spacetime boundary) follows some timelike trajectory 
in a static bulk spacetime. We start off by considering the most general brane-based for- 



malism, following the approach of Ipser and Sikivie ||22|| , finding the general wall solutions 
for a fixed wall of constant spatial curvature embedded in a bulk of constant curvature. We 
then demonstrate how a fixed wall embedded in a non-static spacetime is strictly equivalent 
to a moving wall embedded in a static spacetime thereby establishing in full generality the 
equivalence between the two approaches. This is true for a wall of completely arbitrary equa- 
tion of state separating spacetimes which may even have a different cosmological constant. 
Having shown this equivalence we then proceed to find and study the most general wall 
trajectories, i.e. cosmological evolution equations, in a static AdS or even flat background. 
The reason for studying such a general set-up is that some more recent generalizations of the 
RS scenario have involved not only non Z 2 symmetric walls, such as the 'Millenium Model' 



of Kogan et. al. [2^] consisting of two positive tension walls at orbifold fixed points with 
an additional negative tension wall freely moving inbetween, but also patching together of 
spacetimes with different cosmological constant, such as the GRS model plf , i n which the 
central Z 2 symmetric wall is flanked by two negative tension branes with flat space in the 
exterior. The curious feature of these apparently contrived models is that gravity not only 
changes nature at short but also at ultra-large scales, becoming weaker for the Millenium 
model, and five-dimensional for the GRS model. (In fact gravity for the GRS model has 
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several peculiarities, p5[-[27jj, which may be ameliorated by a hybrid double- wall variant of 
the Millenium model, [25], in which the negative tension wall of the Millenium model is 
replaced by two negative tension walls with a slice of flat space inbetween.) 

While we might expect that cosmology in such models is potentially delicate, the beauty 
of the geometrical approach is that the search for a cosmology becomes a local question 
of finding a suitable trajectory for the spacetime boundaries, thus the presence of extra 
walls is irrelevant - unless these walls collide. We therefore do not consider issues such as 
radion stabilization [p9|, or the like, simply examining the possible trajectories (and hence 
cosmological solutions) of the brane universes. 

In the next section we set up our formalism and then find the general solution to the 
Einstein equations for a constant spatial curvature wall in a constant curvature spacetime, 
generalizing Taub's solutions [[HJ to allow for the presence of a cosmological constant in the 
bulk. We then show how from this brane-based approach we can cross over to the bulk- 
based approach for a wall of arbitrary trajectory evolving in a static spacetime. In section 
III we establish the most general cosmological evolution equations for the different classes 
of spacetime solutions. In the following section we analyze some specific cases of interest in 
cosmology finding the wall trajectories analytically. We make some concluding remarks in 
Section V. 



II. GENERAL WALL SPACETIMES 

In this section we derive the general spacetime of a brane universe, and the equations 
of motion it must obey. As per usual, we shall be modelling our four-dimensional Universe, 
U, to be an infinitesimally thin wall type defect of constant spatial curvature embedded in 
a five-dimensional spacetime. In other words, we are looking for a spacetime with planar 
(or spherical/hyperboloidal) symmetry in three of its spatial directions, which has one (or 
more) hypersurfaces on which the spacetime curvature has a distributional singularity cor- 
responding to a ^-function source for the wall energy-momentum. The most general metric 
admitting this symmetry can be written in the form 

ds 2 = e 2u ^' z \B{t, Z ))- 2 / 3 {dt 2 - dz 2 ) - B 2 / 3 

where k — 0, ±1 represents the spatial curvature of the 3-spatial sections, and u, B will 
satisfy the bulk Einstein equations (with or without a cosmological constant), as well as 
appropriate jump conditions at the wall which we will choose to set at z = 0. 

We start by summarizing the jump conditions for later use. As is conventional, we denote 
the normal to U by n a , in terms of which the first fundamental form, h a b, of the wall is given 
by 

hab = 9ab + n a n b (4) 

which is simply the projection of the bulk metric on U. The second fundamental form, or 
extrinsic curvature of U is 

Kab = h\ a h d b )Vcn d (5) 



dx 2 



1 - KX 



+ x 2 dn 2 n 



(3) 



4 



The intrinsic and extrinsic geometry of IA are related to the bulk curvature via the Gauss 
Codazzi equations 

Racun c n d = K ac K c b - C n K ab = R$> - R cd h c a h d b - K ac K c b + KK ab (6a) 
h bc V b K ca - h b a V b K = R^h b a n c (6b) 
i? (4) - K ab K ab + K 2 = 2G ab n a n b (6c) 

where R^ b is the intrinsic Ricci tensor obtained from h ab , and C n is the Lie derivative with 
respect to n a . 

Since the wall's energy-momentum tensor T™ b is a distributional source, we can now 
integrate (Bal) across the wall, and defining 

S ab = I T2 dl (7) 



we immediately arrive at the Israel junction conditions [31 



AK ab = [S ab ~^Sh ab ], (8) 

where AK ab = K^ b — K~ b is the jump in the extrinsic curvature, and we have set 87rG5 = 1. 

Furthermore taking the sum and difference of (|6cl) and (|6b"D respectively and using (||) 
we have, 

R^ - K ab K ab + K 2 = ^(S ab S ab - ±S 2 ) + 2G nn (9a) 

K ab S ab = 2AG nn (9b) 
V[ 4) 5 6c = (9c) 
<' ] K b a - V^K = (9d) 

where Q = (Q + + Q~)/2 stands for the mean of a quantity across the wall. We should 
stress that (^]) are in fact integrability conditions and hence will result from the Einstein 
and junction conditions (H). 

Since we are looking for cosmological solutions, we shall assume that our brane Universe 
is made of homogeneous and isotropic matter, hence 

S ab = vu a u b + T{h ab - u a u b ) (10) 

where e is the surface energy density and t the tension of the brane U. The timelike vector 
u a is the 5-velocity of an observer comoving with the brane Universe. Obviously if e = t we 
have a domain wall, and if t = 0, a dust wall. 

We begin by finding the most general bulk solution before examining these boundary 
junction conditions. The bulk Einstein equations in this case are simply, 

Rab = -^gah (11) 

where A is the cosmological constant. For the metric @ these give the following system of 
partial differential equations, 
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B 



tt 



B 



zz 



[2AB 1 / 3 - QkB- 1 / 3 ) e 2v 

f,tt - ^,zz = (jj. 

u ,z B ,t + u ,t B ,z = B tz 
2u z B >z + 2v x B x = B tt + fi zz 



13-2/3 + KjB -4/ 3 ^ e 2, 



(12a) 

(12b) 

(12c) 
(12d) 



The particular way of writing the metric (Q) shows the Liouville like character of the system 
of equations, where (|12a|) , ( |12b|) are in the presence of A, non-homogeneous coupled wave 
equations and ( |12c| ) and ( |12d| ) will turn out to be integrability conditions. This is not 
too surprising since we are effectively studying a 1+1 gravity problem by considering the 
Kaluza-Klein reduction of the constant curvature spacelike dimensions. 
It proves easiest to rewrite (|i~2"D in light-cone coordinates, 



u 



in which 



B, 



(2A5 1 / 3 
A 



3-2/3 + KB -m\ e 2u 



2v U B u 



B,u [ln(B,u)] )U 
B >v [ln(B >v )] !V = 2u v B jV 

First note that (|14c| , |14d| ) can be directly integrated to give 

e 2 ^ = V'(v)B >u = U'(u)B >v 



(13) 

(14a) 

(14b) 

(14c) 
(14d) 

(15) 



where V'(v), U'{u) are arbitrary nonzero functions of u, v for generic A, k. (If A = k — 0, 
then it is possible for one of U'(u) or V'(v) to vanish - see below.) It is straightforward 
to see then that B and v have the form (where we remind the reader that 'primes denote 
ordinary differentiation with respect to the unique variable of the function) 



B = B{U{u) + V{v)), 



-2i< 



which reduces the remaining PDE (T^a) to the ODE 

B" - (2AB 1/3 -6kB- 1/3 ) B 



B'U'V 







(16) 



=}► B' - -KB^ + 9kB 2 / 3 = 9fi 

where /1 is an integration constant. This last relation can be integrated to give 

1 



tan - — - — r + coth — 

|r_| ! r + 



V9k 2 - 3//A 

where r = B 1 ^ 3 , c is a constant of integration, and 



(U + V) 



(17) 



(18) 
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Hence the general solution to (JTT]) is given by, 



ds 2 = B'U'V'B- 2 l\dt 2 - dz. 2 ) - B 2 / 3 dx 2 m (20) 



ZH 



where B = B{U + V) satisfies (fP7|) or (|18l), U{u) and V{v) are arbitrary functions and dx 
stands for the three dimensional constant curvature metric written out explicitly in (|]). 
When A = k = 0, ( |TTD gives (wlog) B = U + V , and there are now in fact two allowed 



classes of solution, in the terminology of Taub [p0 |; the above, (0), being a class II solution. 



The class I solutions are distinguished by having U' — or V' — 0, but not both, in which 
case the general metrics are 

ds * = UWg) (dt 2 _ rfz 2 } _ ^2/3^2 (21a) 

rfs 2 = ^^l(dt 2 - dz 2 ) - F 2 / 3 dx 2 (21b) 

for [/' 7^ and V ^ respectively. In fact these class I solutions are flat, as we will see in the 
next section. In summary Einstein's equations (|l|) admit two distinct classes of solutions 
I, ( |21aD , and II, fl2Ti|), which depend on two arbitrary functions. 



Having derived the general bulk solutions, now let us examine the constraints, or bound- 
ary conditions at z = imposed by the wall. For what follows we shall assume Z2 symmetry 
in the z-direction to make notation easier and results more transparent. We shall however 
be dropping this assumption when actually looking for general cosmological solutions in the 
next section. 

Reflection symmetry about z = permits us to consider only z > in (^) since the 
metric is an even function with respect to z. Extrinsic curvature components are thus odd 
functions in z and hence, 

AK ab = 2K+, 7Q b = 

In the coordinate system adapted to the wall the intrinsic and extrinsic geometric quan- 
tities are greatly simplified, and take the form 



E 

— E — t) 



-2e-B~ 2 ^B = -e-» B ' {V '- U ' 




2d z (B^e-) = B^e- " ^) + T^iY* - U')(AB^ - 6„) 



(22a) 
(22b) 



Hence ( p2|) tell us that given two arbitrary functions U and V we can determine e and t. 
We now show that only one of these functions is physical. 

As a general rule any coordinate transformation in the bulk will shift our boundary 
making our wall evolving in time. However we can make a class of coordinate transformations 
in the bulk which leave our wall boundary conditions fixed, namely 

u — /(it), v — f(v) (23) 



i.e. the rescaling of the light-cone coordinates by some arbitrary function. This boundary 
invariance is a result of the two-dimensional conformal symmetry on the t-z plane. Now 
we choose a particular coordinate transformation, / = V, to fix this non-physical gauge 
freedom, thereby effectively setting V = v. Hence there is only one remaining arbitrary 
function or physical degree of freedom, U(u), which completely determines e and t for the 
wall. 

There are now two equivalent ways, as it turns out, to look for cosmological solutions. 
We either, in the spirit of [|T2j| , consider a fixed boundary in the general spacetime (p0|), 
using the junction conditions fl2"2"|), or we make coordinate transformations simplifying the 
bulk and study the moving wall trajectories in that case. 

In fact, it turns out that the general class II solution ([H]) is a static spacetime. Setting 

r = B 1/3 ; t = 3(v~U) (24) 

(^0|) becomes, 

dv 2 

ds 2 = h(r)dt 2 - — — - r 2 dx 2 JU (25) 
h(r) 

where h(r) has the familiar form 

B'B- 2 ^ A „ n 



Mr) = — k — — v — —r (26) 

V > 9 6 r 2 y ' 

recognised as the generalised "Schwarzschild-AdS" solutions. 

Note however that in making our background static our wall is no longer fixed, since 
the trajectory z = 0, or u = v , gives a nontrivial relationship between r and t: r — R(t) 
from (|24]) . The reader should note here, that time dependence of the metric fl2"0|) in the 
bulk, portrayed by arbitrary function U(u), and static boundary (^) has been transformed 
into time dependence for the boundary portrayed by the arbitrary trajectory r = R(t), in a 



static bulk spacetime (p5|) . Hence all essential information of the bulk metric is transferred 
onto the boundary conditions and vice-versa. 

As a final remark, we note that the generalization to domain walls in arbitrary dimension 
is straightforward, simply altering the fx/r 2 to yu/r^ n ~ 3 ^ in ( PB"|) (see JJIJ for a different 
approach) . 



III. COSMOLOGICAL EVOLUTION EQUATIONS 

In this section we derive the equations of motion for a general (i.e. not necessarily sym- 
metric) wall bounding two regions of spacetime with possibly differing cosmological con- 
stants. The reason for setting up the problem in its full generality is that we will be able 
to use the same set of equations to deal not only with non Z 2 symmetric walls, such as 
in the Millenium model, but also with walls which separate negative and zero cosmological 
constant spacetimes. This latter case will of course also require consideration of both class 
I and class II solutions for the A = k = spacetime. We begin with a wall bounding two 
class II spacetimes. 



S 



The starting point is the general solution (p5|). To get a wall solution, we simply take a 
general boundary 



X a = (t(r),R(r),x,e,4>) (27) 

and compute the extrinsic curvatures on each side. Note that r is the proper time with 
respect to an observer comoving with the wall, i.e. 

i 2 h(r) -R 2 h~ l (r) = 1 (28) 

and is therefore well defined and intrinsic to the wall, whereas t is the coordinate time in 
the bulk, and need not agree with the coordinate time on the other side of the wall. The 
only other subtlety arises in the choice of sign for the outward normal 



71 n 



±(R(r),-i(r),0) (29) 



This corresponds to the fact that we have a choice of which part of the spacetime we wish 
to keep (i.e. r < R or r > R). This depends on the physical solution required, i.e. whether 
we have a positive or negative energy wall, and what sort of asymptotic solution we require. 
In general, positive energy walls match across two interior spacetimes and correspond to the 
lower (— ) choice of sign in the outward normal and vice versa. It is also possible to have a 
wall matching an interior and exterior spacetime (for example the spherical collapsing wall, 



22| , used by Gogberashvili |33|] in an early version of the RS scenario), however, for clarity 
in what follows unless explicitly stated, the formulae will be valid for +(— ) energy walls 
matching two interior (exterior) patches of the Schwarzschild-ADS spacetime. 
Computing the extrinsic curvature components on each side of the wall gives 

E = T^(i+h++Lh-) (30a) 

where t± represent the coordinate times on each side of the wall, i±h± being given by 
(f28|) in each case. Note that we can also express ( |30b|) in a different fashion as an energy 
conservation equation, 

De R , , 
- + 3-( E - T ) = (31) 

These represent the general equations of motion for an arbitrary wall separating two 
regions of spacetime with possibly different cosmological constants. Note that at this point 
e and t are unconstrained physically, and are simply given by the above expressions in terms 
of -R(t), which, other than representing a non-spacelike hypersurface, is completely arbitrary. 
This is worth stressing - that we can take any such R(t), and we will get a wall solution with 
energy and tension determined by the embedding. Only when we impose physical constraints 
on our energy and tension do we get constraints on our possible hypersurface R{t). Before 
exploring such constraints however, we simply note that (|30|) can also be expressed in terms 
of brane cosmological evolution equations (cf. Binetruy et. al. [|12|]) which are 
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R\ _ e 2 + 9(Afc 



R 



36 



t 9Afc 2 A/i 9(A/i) 2 



(32a) 



R= _(2e-3t) 



f + gf^g) + R f ' ^) _ f) - A _ g^M! (32b) 

36 2E 3 i? 3 / \ 4e 3 I R 3 4E 3 i? 7 



where we have now put 6&; 2 _ = — A±. These are the completely general brane cosmological 
evolution equations for a brane separating two 'ADS' regions, and were first derived by 
Stoica et. al. [|4[ . 

Since we are particularly interested in the cosmological solutions for branes in models 
with regions of vacuum spacetime, we must also consider the case of a planar wall bordering 
on a class I spacetime which we may write as 



ds 2 = H{y)drdv — r 2 dx 2 



(33) 



This spacetime is actually flat [30], as can be seen by the coordinate transformation 



x = rx 

t* -z* = r 

t* + z* = [ Hdv + rx 2 



(34a) 
(34b) 

(34c) 



which gives the Minkowski metric in the starred coordinate system. 

Clearly to match with an ADS spacetime we need to set r = R(r), and compute the 
relevant extrinsic curvature quantities. For a wall bounding a class II and class I spacetime, 
we find that the junction conditions are now 



R 



R + i-h-) 



— E 

3 



± (R R+W- 
\R Lh- 



(35a) 
(35b) 



Therefore we now find a somewhat different set of cosmological equations for the evolution 
of the wall: 



R 
R 


\ 2e 




2> 2 - 



E 3/x_ 
6 ~ 2e,R 4 



(2e - 3t)R 
36b 3 ' 



81k A - 



li- 



2e 3 R 3 



9k 2 (3r - 2e) - e 3 



27t/x 2 

4E 3 i? 7 



(36a) 
(36b) 



where we see that the 'Friedman' equation is now linear in R. 

Finally, if we are matching two class I spacetimes, the junction conditions give the 
particularly simple relations 



6R 

( AR 2R' 



(37a) 
(37b) 
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In the next section we will compute cosmological trajectories for general walls, however, 
before doing so it is useful to verify that this approach is indeed valid by cross-checking 
it against a few simple known solutions. For example, a planar symmetric domain wall in 
pure Einstein gravity is known to have a particularly simple interpretation in terms of the 
matched interiors of two hyperboloids in Minkowski spacetime ||35|| . In the context of these 
equations, a A = k = spacetime can match class I or class II spacetimes. For the matching 
of two class I spacetimes we must use ( |3"7| ) with e = t. Clearly this has the solution R = e Er / 6 , 
giving a de-Sitter like induced metric in agreement with the standard domain wall metric 
in wall-based coordinates (see e.g. P2"| , |3"6]] ). Inverting this relation to find R{v) in terms 
of the bulk coordinates gives R = —18/e 2 v. Finally, inverting ( |34"D gives the particularly 
simple bulk equation of motion for the domain wall trajectory: t* 2 — z* 2 — x* 2 = — 36/e 2 , i.e. 
a hyperboloid in Minkowski spacetime. The 'planar' domain wall is therefore a boundary 
between the interior of two Minkowski hyperboloids as required. 



IV. COSMOLOGICAL WALL SOLUTIONS 

In this section we apply the general equations worked out in the previous section to 
find a variety of cosmological solutions. Let us begin by considering the particularly simple 
equation of state of a domain wall, e = t, and find the most general trajectories. First of all, 
note that (|31~1) implies that e is a constant, therefore the brane-Friedmann equation, fl32a|) . 
can be written as 

D4 / A T3l \ 2 

= aR 8 - kR 6 + bR 4 + c (38) 




where 



+ ^-F (39a) 



36 4 



b = T2 (39b) 

(39c) 

4e z 

which can be integrated in general, giving R 2 in terms of elliptic functions, although the 
expression is not particularly illuminating. If one of k, c is zero, the solutions are very simple 
to write down, but for general a, b, c, k, the qualitative behaviour of the solutions can be 
deduced from (j38|) . The constant 'a' can be seen to measure the departure from criticality 
of the domain wall - i.e. that value which allows a static planar domain wall solution in 
the absence of an ADM mass term. We will therefore denote a = walls as critical, and 
a > (<) walls as super-(sub-)critical. 

Clearly, if there is a nonzero ADM mass in the bulk there is always a solution which 
expands outward from an initial singularity, R — 0, which is where the wall trajectory 
touches the central singularity of the black hole. Whether or not this solution expands 
indefinitely, or there is a final singularity depends on the roots of (|38|). For a > 0, there are 
zero or two (possibly repeated) roots for R 2 > 0. If b > jj- for k — 1, or b > otherwise, 
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then there are no positive roots and the cosmology expands indefinitely. Otherwise there can 
be two separate roots, in which case there is a cosmology with initial and final singularities, 
and a nonsingular cosmology which contracts in to a minimum radius and re-expands. If 
there is a repeated root, i? 2 , then the cosmology asymptotes R c exponentially at late times 
either as a contracting branch, or an expanding branch with an initial singularity. If a < 
there is always one root and our universe starts and ends its life on the black hole singularity. 
Finally, there is always a static solution at any positive root of the quartic. 
For planar domain walls, i.e. k = 0, we have the exact solution 

JL (g^r-To) _ 2b + {b 2 - 4ac)e- 4 ^( T - T °)) a > 

RA =\ %^n{4v^(T-T )}-2H a <° (40) 

4b(r - r ) 2 - f a = 

which clearly demonstrates the expansion outwards from an initial singularity, and for a > 
shows the late time inflationary nature of the cosmology. 

If c = 0, i.e. the mass terms on either side of the wall are the same, the exact solutions 



arc 



_L (g2VS(T-T ) + 2k+ (k 2 - Aab)e- 2 ^ r - T ^) a > 

R2= \ ^ Msin { 2 V^( r - r o)} + ^ a<0 (41) 

-k(t - r ) 2 + £ a = 



where we have taken k ^ 0: if k = use (|40|). A nice analysis of supercritical domain wall 
solutions and their bulk embeddings was given by Khoury et. al. in [j37H - 

For the case of a general wall we need some sort of equation of state t(e), for which we 
will follow convention in setting 

e = e + p ; T = E -p (42) 

where p = (7 — l)p would be a typical equation of state, although we re-emphasize that 
for p — e there is no reason to suppose that such an equation of state will be accurate, or 
even possible. The conservation of energy momentum ([31]) implies that p = po-R -37 and 
hence e is no longer a constant. In order to deal with this additional difficulty (in the 
non-Z2-symmetric case) we keep only terms linear in p which boils down to looking at wall 
trajectories for late times. For example in the case of radiation cosmology, 7 = 4/3, the 
Friedmann equation (|32ap reduces to, 



R\ _ e 2 9(Afc 2 ) 2 - k 1 /_ 9AA: 2 Ap Eo p 9p (Afc 2 ) 2 
R " 36 4e 2 R 2 + R*\ fl 2e 2 + 18 2e§ 



(43) 



The above equation has exactly the same solution as PH ) for k 7^ where however the 
constants a and b are given by, 

„ _ 4 1 9(AA; 2 ) 2 F . 9APAp e oPo 9 Po (AA: 2 ) 2 
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Another important class of solutions are those with symmetry, for which the cosmo- 
logical Friedman equation has the form of (|38|) with the differences 'A' of all quantities being 
zero. The Friedmann equation ( |32a| ) becomes 

^V=— -A- 2 - — 4- EoPo 4- JL + Po (A^) 
\R J 36 R 2 18R 3 ~< i? 4 1 } 

where the above is now an exact expression. For example, a radiation cosmology would have 
7 = 4/3, and the equation of motion would be as (ffSD, but with the constants a r ,b r ,c r , now 
defined as 

2 2 
E ; 2 l , E 0P0 Po /i fi \ 

a r = k ; o r = u H ; c r = — 46 

36 P 18 ' 36 v ; 



and so a generic planar solution would be given by ( f40D with these new values for the 
constants a, b, and c. We see now the interplay between the contribution of the cosmological 
matter energy density p , which contributes to both b and c, and the "CFT" contribution, 
/i, which appears only in b. This changes the way in which the initial singularity, R = 0, 



is approached, for example, Gubser's solution has a = = p , and R = /i*y2(r — r ), as 

opposed to R = fe3 [4(r — r ) 2 — c/b 2 ]^. A description of Z 2 symmetric cosmological solutions 
and their embeddings was given by Mukohyama et. al. WJfi . 

Let us now turn to the less well-explored question of cosmological solutions in models 
with slices of vacuum and ADS spacetimes, such as the noncompact GRS model, or compact 



models such as those in p8| , |38|] . Without loss of generality, we will take k + = 0, = k, 
for which the general equations of motion are given by or (|3"6|). For simplicity, since 
we are unlikely to be interested in cosmologies on negative tension walls, let us consider 
the case where the boundary has the equation of state of a domain wall, i.e. e = t, and 
let e = — (e + 3k) be the departure from 'criticality' of this domain wall (e > being 
supercritical). 

We start by considering the matching of two class II spacetimes. Note that for k = 
0,-1, in order to have the correct spacetime signature, we require /i + < 0, and we have a 
naked singularity in the full bulk solution. For k = 1, p + > gives the five- dimensional 
Schwarzschild solution. 



Substituting the values of k± gives for the constants in 

e 2 , , 2 £(£ + 2e)Au 9(Au) 2 . . 

Note that unlike the generic domain wall case, a > for both sub and super-critical walls. 

For k = 1, > and we again have solutions which expand out (or contract into) 
the singularity at r = 0. Depending on the roots of the quartic, there are solutions which 
expand indefinitely, recontract, or indeed asymptote a constant radius. The critical domain 
wall always recollapses. For k = —1, /i + < 0, and the qualitative behaviour is similar, 
however the critical domain wall now expands indefinitely if 27c > — 4/x + . 

For the planar (k = 0) walls, the exact solutions are 

^4 _ / Ta (e ±4 ^( T ~ T 0) - 2b + {b 2 - 4ac)e T4 ^ T - T0 ^ a > ( . 

7= \4b(T-T ) 2 -l a = [ ' 
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FIG. 1. A negative tension domain wall patching vacuum to Schwarzschild-ADS spacetime. 
The vacuum bulk has a naked singularity. 



For the critical domain wall, b = /i + < 0, and we see that the domain wall universe has 



an initial and final singularity, for r — jj^- < r < r + jjp-- Inverting (|28|) to find the 



trajectory in terms of coordinate time in the vacuum spacetime to the right of the wall shows 
that these singularities occur at finite coordinate time, and the wall actually expands out of, 
and re-collapses into, the central naked singularity in this spacetime. This corresponds to 
figure HJ, and would appear to be a rather undesirable spacetime from the five- dimensional 
point of view. 

For the subcritical domain wall however, b < with b 2 — 4ac > 0, hence the cosmology 
has no initial singularity, but merely follows a modestly corrected sinh trajectory. For the 
supercritical domain wall, there are two possibilities, depending on the relative magnitudes 
of /i±. If 9k 2 \fi + \ > fi-(e 2 + 6ke), then the cosmology is completely nonsingular as for 
the subcritical case, however, if 9k 2 \fi + \ < fi-(e 2 + 6ke), then there is an initial (or final) 
singularity with the universe inflating away (deflating towards) it. 

Finally, for a planar asymptotically vacuum spacetime, the other possibility is matching 
to a class I spacetime for which we need fl36|). Setting e = t = — e — 3k < again we have 

(49) 



R e(e + 2e) 



R 6e 2e# 4 

which has the general solution 




R * = , " ^ I 3E V < e(e+2E) ~ T » (gQ) 

- t ) e = 

A supercritical wall therefore has R > and expands outward from an initial singularity 
(where it touches the Schwarzschild singularity on the (— ) side of the wall). A subcritical 
wall however has R^O for -R 4 ^ r-r^r - This shows that there are two solutions, one of which 

< > e(e+2E) ' 

has an inital singularity and expands outward to i? 4 = £ ^2e) > anc ^ an °ther which contracts 
to R c . A critical domain wall has a power law expansion outward from an initial singularity. 

It is interesting to transform into the starred coordinate system to follow these wall 
trajectories in Minkowski spacetime: 
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FIG. 2. A negative tension domain wall patching a class I (flat) vacuum spacetime to 
Schwarzschild-ADS. The dotted line indicates a possible additional positive energy cosmological 
domain wall. 



x* = i?(r)x (51a) 
= R( T ) (51b) 

t* + z* = I + i?(r)x 2 = F{R{t)) + i?(r)x 2 (51c) 
J R{t) 

For a critical domain wall F(R) = 4k 2 R 7 /7 '/z 2 , and for a noncritical domain wall bordering 
pure ADS, F(R) = —1/4R. We see therefore that this latter case is the familiar hyperboloid 
in Minkowski spacetime, although it is now the interior which is excised for this negative 
tension wall. 

The critical domain wall satisfies 

Ah 2 

x *2 = f 2 _ z *2 _ _ ^8 (52) 

which is a deformed hyperboloid. At fixed time, the 3-sphere is now squashed into an egg- 
shape in the |x*|, z* directions. Again, it is the interior of this squashed S 3 which is excised. 
On the ADS side, the wall expands outward from the initial black hole singularity and at 
late times expands as R oc t_ (c.f. R oc t for Gubser's critical Z 2 symmetric solution). The 
appropriate matching is shown in figure 

Finally it is worth stressing that these trajectories are independent of whether or not 
there is an additional positive tension wall present - the existence of the wall trajectory is a 
local question. We can therefore add in a positive domain wall in the usual way by simply 
adding another boundary to the Schwarzschild ADS spacetime as shown by the dotted line 
in figure ||. 



V. DISCUSSION 

Starting from the most general five- dimensional metric with homogeneous and isotropic 
spatial 3-sections, we have shown that the most general cosmological brane solution matches 
two Schwarzschild-ADS spacetimes, or a Schwarzschild-ADS spacetime with a class I or II 
vacuum spacetime. It is worth emphasizing that the three dimensional homogeneity and 
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isotropy ensures that the bulk metric solution is invariant under two dimensional conformal 
symmetries. By inputting a boundary or wall we in principle break these conformal symme- 
tries, however, it turns out that the junction conditions are such that only half are broken. 
This in turn ensures that the only remaining degree of freedom is the wall's trajectory itself. 
Hence all the essential information of the wall's dynamics is contained on the wall trajectory 
itself (the boundary) the bulk being a specific static spacetime playing a 'background' role. 

We then derived the generalised FRW equations for a wall with general energy and 
tension, then applied these to a range of cases, focusing on the previously unexplored case 
of a wall bordering a vacuum region on one side. The interesting feature of these latter 
solutions is that if two class II spacetimes are matched, the bulk can contain timelike naked 
singularities on the vacuum side of the wall. For a more satisfactory nonsingular solution, 
one must match the planar wall to a class I spacetime. Here the wall trajectories are 
reminiscent of the vacuum domain wall of Vilenkin, Ipser and Sikivie, in that they are 
deformed hyperboloids, although it is the interior of the hyperboloid which is excised for 
this negative energy wall. 

If we now wish to construct a 'cosmological' GRS solution, i.e. one which has a Z2 
symmetric positive energy central wall with matter residing on it, and a negative energy 
outer wall, then we must combine one of the Z 2 solutions of (H5I) with a class II /II or class 
II /I solution for the outer wall. In each case, we must be careful to keep the (— ) wall to 
the right of the (+) wall, R_ < R + . For example, if we take the planar domain wall, the 
central domain wall trajectory, R + , is given by fl40|) with the constants a, b, c in ([IB]). Using 
the class II/I solutions of ( 48|) for the (— ) wall shows that we cannot match arbitrary walls, 



since for example a subcritical central wall would collide with a critical or supercritical outer 
wall bringing our universe to an abrupt and catastrophic end. Two generic critical walls are 
compatible however, with 




36 M +12fcpo (53) 

as is a critical central wall with a subcritical outer one. A supercritical central wall can exist 



with a critical or subcritical outer wall, or a supercritical outer wall with 2y'E^_ — 36fc 2 > 
(9k 2 — e 2 )/|e_|; and a subcritical central wall can exist with a subcritical outer wall provided 

yjbl-4a r c r -b r 9^ 



^ {9k 2 -E 2 _)- 

For example we can match the critical Gubser solution, R^_ = 4^r?, with the critical 
class II/I solution R*_ = 2fir_/k. However, a pure critical planar radiation cosmology, 
R+ = p (l6k 2 T+ — l)/12k cannot be matched to any class I asymptotically flat spacetime 
across a (— ) wall, as the only possible solution is an exponentially expanding or contracting 
one for a super /sub critical wall respectively. 

Therefore, while there are some constraints on the central wall cosmologies, these are 
not overwhelmingly restrictive. 

This leads us to the issue of the radion in these cosmological models. There has been a 
great deal of debate about the radion in the context of cosmological models (see e.g. p9| ) 
and in the case of quasi- localised gravity (see e.g. [p5H27fl). We have not supposed any sta- 



bilization mechanism, [B9 , for the extra dimension, but simply looked for free cosmological 
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solutions. Since the radion leads to anti-gravity in the original GRS model, it is interesting 
to see how it manifests itself here. 

and it 



The wavefunction of the radion is straightforward to find in linearized gravity, [40 



behaves as a scalar field living on the brane worldvolume which satisfies a massless equation 
of motion, d 2 f = 0. In the context of a homogeneous and isotropic cosmology, this has the 
simple solution / = fot, where t is the coordinate time on the brane, and we have chosen to 
have zero displacement for t = 0. Now consider a supercritical ^-symmetric positive energy 
planar domain wall in pure ADS spacetime; this follows a trajectory 



R = expjvV - 36k 2 (r - r )/6} 



\ k Ve 2 - 36Jfe 2 , 

1 (t-t ) 



(54) 



in terms of either the brane proper time, r, or the local bulk coordinate time, t. We see 
therefore that for small t and Ve 2 — 36k 2 , 



1 + 



kyZkyiz — 6k 



(< - to) 



(55) 



we can therefore identify /o = & 3//2 \/e — 6ky/3, i.e. 5e = e — 6k = 3/qA;~ 3 / 2 , which has the 
correct dependence on / for a standard energy momentum tensor for the radion. (Note 
we cannot have a subcritical planar domain wall.) This means that a negative energy wall 
requires e < —6k, and nominally a negative radion energy on that wall. 

Obviously these solutions do not see any evidence of the fifth dimension "opening up" 
p4| , or any instability due to the radion, however, this is because of the high degree of 
symmetry of the set-up. If we were to break isotropy or homogeneity, the integrable nature 
of the system would likely be destroyed, and the simplicity of the description of the wall as 
a trajectory in bulk Schwarzschild-ADS would disappear. Destruction of these symmetries 
would lead to a brane-bulk interaction (as evidenced by the extreme case of the as yet 
undiscovered black hole on the brane solution - see for discussions on this topic) through 
which the more familiar effects of the radion would be recovered. 

Of course it is tempting to enquire whether the so-called "missing mode" of free motion 
of a wall can be identified in this set-up. For a non-gravitating defect, the effective equation 
of motion is given by the relation K = 0, i.e. the wall (in this case) is a minimal hypersurface. 
An equivalent perturbative description would state that the displacement of the wall satisfies 
a massless wave equation on the wall. There have been claims El in the literature that 
when gravity is included this 'free motion' disappears, as would appear to agree with the 
rejection of the free 'radion' solution 



4k 



-K 



[IV 



(56) 



in the usual Randall- Sundrum perturbation analysis. However, in this simple picture, each 
wall has its own 'radion' corresponding to motion through the ambient spacetime. Indeed, 
a perturbative analysis of the thick four- dimensional domain wall [43] indicates no such 
disappearance of this motion with the coupling to gravity. What therefore can be going 
on? A clue perhaps lies in the effective equation of state for a freely moving wall. As 



Carter and Vilenkin ||44|| , argued for the cosmic string, a freely moving defect will have an 
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altered effective equation of state. The effective energy per unit area will increase, and the 
tension will decrease. Generalizing the Carter- Vilenkin formula for the domain wall gives 
the effective equation of state 



,3 



ET 



(57) 



for a 3-brane domain wall. Clearly for e close to e , this gives an equation of state of a 
radiation cosmology on the background domain wall. Now we see how a freely moving wall 
might no longer be a small perturbation of a 'straight' one: the averaged equation of state of 
the freely moving wall causes the wall to follow a centre of mass trajectory associated with 
a radiation universe, which is not, at late times (when the radiation universe approximation 
is particularly good) a small perturbation of the unmoving wall. Of course, a freely moving 
wall will not be isotropic and homogeneous except at the very large scale, and it is also 
possible that there are brane-bulk interactions which further complicate the issue. 

Finally, the problem of brane cosmological perturbation theory is very important to 
understand if we wish to do real cosmology. The setting up of the formalism for perturbation 
theory on the brane is already underway |]45|j46[| ; by considering the cosmological domain 



wall in terms of its global spacetime structure, the role of the bulk, its interactions with 
matter on the brane, and the interpretation of some of the gauge invariant variables in 
terms of bulk physical quantities will be more clearly elucidated. 
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